A nonlinear scalar model of extreme mass ratio inspirals in effective field theory 

I. Self force through third order 

Chad R. Galley* 1 

^ Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, 91109 
'^Center for Fundamental Physics, Department of Physics and 
Center for Scientific Computation and Mathematical Modeling, 
University of Maryland, College Park, MD 20742 

The motion of a smah compact object in a background spacetime is investigated in the context 
of a model nonUnear scalar field theory. This model is constructed to have a perturbative structure 
analogous to the General Relativistic description of extreme mass ratio inspirals (EMRIs). We apply 
the effective field theory approach to this model and calculate the finite part of the self force on 
the small compact object through third order in the ratio of the size of the compact object to the 
curvature scale of the background (e.g., black hole) spacetime. We use well-known renormalization 
methods and demonstrate the consistency of the formalism in rendering the self force finite at 
higher orders within a point particle prescription for the small compact object. This nonlinear 
scalar model should be useful for studying various aspects of higher-order self force effects in EMRIs 
but within a comparatively simpler context than the full gravitational case. These aspects include 
developing practical schemes for higher order self force numerical computations, quantifying the 
effects of transient resonances on EMRI waveforms and accurately modeling the small compact 
object's motion for precise determinations of the parameters of detected EMRI sources. 

I. INTRODUCTION 

The Laser Interferometer Space Antenna (LISA) [Tj is expected to see thousands of low-frequency gravitational 
wave sources when it is launched around 2025. LISA's detection of extreme mass ratio inspirals (EMRIs) - a white 
dwarf, neutron star or small black hole (collectively referred to as small compact object, SCO) inspiraling toward a 
supermassive black hole of mass 10^ — 10^ solar masses - is expected to provide an unprecedented level of insight into 
the structure of spinning black hole spacetimes as well as into the dynamics and populations of EMRIs in galactic 
nuclei. 

EMRI signals will be extracted with matched filtering techniques from LISA's data stream. Once the presence of 
a signal is established one can use the resulting roughly-estimated parameters to refine the search using waveforms 
calculated with accuracies better than one cycle in the roughly 10^ cycles that are expected to accumulate during 
the last year of inspiral j2j . These high-accuracy templates are expected to be sufficiently accurate that the masses 
can be determined to about one part in 10'* [3 and will require modeling the motion of the SCO with fractional 
accuracy of 10~^ or better. The SCO's motion and corresponding gravitational wave emission may be described using 
a perturbative treatment since the mass ratio for EMRIs is very small, between about 10"^ and 10~^ for EMRIs in 
LISA'S detectable bandwidth. 

Gravitational waves from EMRIs carry energy, linear momentum and angular momentum that results in a force 
on the compact object, called the self force, which causes the SCO to inspiral toward the supermassive black hole. 
The self force on the SCO is conceptually different than, for example, the radiation reaction on a point charge in 
flat spacetime. Radiation reaction is local in time and is typically proportional to the time-derivative of the charge's 
acceleration. However, the self force also accounts for the history-dependent force arising from the interactions of the 
SCO with waves emitted in the past that have backscattered off the background spacetime curvature. As such, self 
force is nonlocal in time and depends on the past motion of the SCO. The self force can be computed perturbatively in 
powers of the mass ratio or, more generally, in powers of the size of the SCO to the curvature scale of the background 
spacetime it moves on. To develop templates with sufficient accuracy to achieve LISA's science goals necessitates 
corrections through at least second order in the mass ratio. We motivate this last statement with a more detailed 
discussion below. 

In this paper, we take first steps towards building high-accuracy EMRI waveforms by deriving the expressions for 
the self force on a SCO due to the emission of scalar perturbations. Specifically, we calculate the self force through 
third order in the mass ratio within a class of nonlinear scalar models that is constructed to be analogous to the 
kinematical structure of the perturbative General Relativistic description of EMRIs. 
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A. High-accuracy EMRI waveforms 

Current theoretical techniques for generating high-accuracy waveforms are significantly underdeveloped when com- 
pared with those methods producing less accurate waveform templates, which are more useful for detecting EMRIs 
(for a review see [HIS]). The former waveforms include "Capra" waveforms (see e.g., and "two-timescale" waveforms 
O [7] . Capra waveforms constitute the highest standard of accurate source modeling for EMRIs since the underlying 
calculations are based on a minimum number of assumptions and are thought to represent the binary's evolution and 
gravitational wave emission most accurately. These waveforms are sourced by the solutions to the self force equations 
of motion describing the perturbed motion of the SCO in the background supermassive black hole spacetime. The self 
force equation through first-order in the (very small) mass ratio is a complicated integro-differential equation, called 
the MiSaTaQuWa equation [8(i9j, for the SCO's worldline coordinates. Accordingly, not a single Capra waveform has 
been computed despite recent advances and progress with numerical computations from several research groups; see 
Ref. [TU] for a recent review. 

It was argued in 112) that if only the first-order (dissipative part of the) self force drives the quasi-circular 
inspiral of a SCO of mass m moving in a Schwarzschild spacetime with mass M, then the accumulated phase of the 
gravitational waveform over an inspiral time ~ M/e, where e = m/M <C 1, is, schematically, 

$-^+0(e°). (1.1) 
e 

The last term represents the error in the phase from not including second and higher-order self force effects, which 
represents an 0(1) correction. Thus, to produce waveforms accurate to less than a cycle requires that second-order 
self force corrections be included in determining the SCO's motion. 

Two-timescale waveforms are based on a systematic adiabatic expansion in which the typical orbital period Toib is 
small compared to the inspiral timescale Tjnsp. To leading order in Torb/Tlnsp, Mino |6j showed that the waveform 
phase depends only on the time average of the dissipative part of the first-order self force. Hinderer and Flanagan |7] 
extended Mino's work by placing it within a systematic two-timescale expansion to calculate "post-adiabatic" (PA) 
corrections to Mino's result. They found at IPA that the time averaged dissipative part of the second-order self force 
is just as important as fiuctuations in the conservative part of the first-order self force. Therefore, second-order self 
force corrections are important to maintain the consistency of the inspiral's adiabatic evolution. 

B. Transient resonances 

Recently, Flanagan and Hinderer [T2] discovered that the SCO may undergo transient resonances during the course 
of its inspiral. Resonances occur only when the SCO evolves on a non-equatorial and eccentric orbit in a spinning 
supermassive black hole spacetime. Such orbital configurations are expected to be generic for EMRI sources detectable 
by LISA. During a transient resonance, the frequency of the true gravitational wave signal undergoes sudden jumps and 
an adiabatic treatment of the inspiral breaks down. Thus, the corresponding signal-to-noise ratio may be significantly 
diminished if using a template bank of two-timescale waveforms. 

One way to address this problem is to provide increasingly accurate waveforms that are capable of tracking the phase 
evolution of the SCO even through the transient resonances. There are at least two ways to do this. One may patch a 
kludge waveform [HI [T5] (describing the resonant phases) into a very accurate two-timescale waveform (describing the 
system at all other times, during the adiabatic inspiral phases) or one may develop Capra waveforms that incorporate 
higher-order self force corrections to the SCO's motion. The former approach is a reasonable possibility since the 
durations of the transient resonances are short compared to the radiation reaction time scale but will be only as 
accurate as the kludge model used. The latter method is the most direct and accurate way but may be difficult to 
realize given that even first order Capra inspiral waveforms have not yet been computed. Nevertheless, in the presence 
of resonances the accumulated phase of the waveform is, schematically, 

$-- + ^ + 0(£°) (1.2) 
e y/e 

and the 0{e^^^^) contribution, which originates from passing through a resonance, requires knowing a part of the 
self force at second order in e |T3]. Therefore, for EMRIs that pass through at least one resonance, the second and 
possibly higher-order contributions to the self force are especially important. 

Estimates of the effects on the inspiral waveform phase from passing through a resonance indicate that the phase 
may change by ~ 20 rad for an EMRI with a mass ratio of 10~^ [13_. Not only is this a significant change for one 
resonant crossing but the effect accumulates for each resonance encountered during the inspiral. As a result, even 
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detecting EMRIs with LISA could be affected by the dephasing from transient resonances. Thus, the motion of the 
SCO is needed with (possibly very) high accuracy in order to sufficiently describe its evolution before, during and 
after each transient resonance. 

Based on these previous works, there are several indications suggesting a need to model EMRI sources with high 
accuracy, which will require incorporating the effects on the SCO from second (and possibly higher) order self force 
corrections. 

C. Intermediate mass ratio inspirals 

Higher order self force corrections may also be needed to model binaries with less extreme mass ratios and could 
be useful even for those with comparable mass ratios. It is natural to think that including higher-order self force 
corrections in the SCO equations of motion will allow for the mass ratio to be relaxed to higher values. Doing so may 
offer the only way of describing binaries with intermediate mass ratios (IMRs) , which have mass ratios in the range of 
~ 10""'^ — 10~^, since neither the post-Newtonian approximation nor numerical relativity are particularly good tools 
for studying the inspirals of IMRs (however, see the recent work of [16]). In fact, an alternative approach based on 
self force methods for binaries with IMRs may be useful for calibrating semi-analytical models (e.g., Effective One 
Body fTT]) and phenomenological hybrid waveforms |18| . 

D. Self force in scalar models 

In this paper, we introduce a class of nonlinear scalar theories that is developed with a structure very similar to 
the perturbative General Relativistic description of EMRIs and will serve as a scalar analog of these sources. We 
calculate the finite (or regular) part of the scalar self force on the SCO through third order in the ratio of the size of 
the SCO - Rm - to the background curvature length scale - TZ - and denoted hy e = Rm/Ti-, which is just the mass 
ratio in the strong-field regime of a supermassive black hole. We take the background to be specified for all time and, 
for simplicity, do not include the effects from the scalar field's stress-energy on the spacetime so that the scalar field is 
analogous to the propagation of metric perturbations on a fixed background spacetime. In regularizing the formally 
divergent self force expressions we use the standard tools of renormalization borrowed from the fields of high energy 
physics and condensed matter. The methods used in this paper will be of direct use for calculating the second order 
gravitational self force for EMRIs in a future paper. 

Historically, scalar models offer a simpler framework for studying the underlying issues of self force regularization 
and for developing practical self force computational schemes. Indeed, the most useful regularization scheme (both in 
terms of physical intuition and practical computations) was first developed and understood in the context of the self 
force on a scalar charge from a linear scalar field in [12]. The first numerical computation of the self force (which was 
performed for a circular geodesic in Schwarzschild spacetime) was accomplished in a linear scalar theory in |20| and 
predated, by about seven years, the corresponding computation in the gravitational case |21| . 

Because of the relative simplicities that scalar models afford, it seems likely that the physics of higher-order self 
force effects can be investigated more easily and quickly than in the gravitational EMRI context. Specifically, one 
can address the qualitative and quantitative effect that higher-order self force corrections have on the waveforms 
themselves (i.e., the change in the phase compared to first-order accurate waveforms, parameter estimation, etc.); 
how transient resonances and the number of resonances encountered during the inspiral affect the waveform and the 
SCO's motion; and how much the mass ratio can be relaxed to higher values while still having a reasonably accurate 
description of the system (the accuracy can be addressed using the next higher order to bound or estimate the errors) . 
It is likely that these questions can begin to be studied using the results from this paper and with small modifications 
to some self force codes currently in use (particularly those using 3-1-1 methods as in [22,^ since these do not rely on 
a mode decomposition of the field) . 

E. Previous work in higher-order self force corrections 

The first work (that we are aware of) regarding higher-order self force computations was carried out by Burko in 
Ref. [lT\. He computed the second-order self force on a scalar charged particle in quasi-circular orbit in a Schwarzschild 
background. Burko derived a formal expression for the equations of motion on the particle for this scenario. The force 
on the particle included contributions from the second-order self force and also from the product of two first-order 
pieces. Since the second-order expressions were unknown, Burko retained only the latter contributions and estimated 
the change in the accumulated phase due to these self force corrections when compared to the first-order accurate 
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waveform phase. He found that there was a relevant correction of 0(1) cycles to the first-order accurate phase. The 
change to the phase from amplitude corrections was found to be about a tenth of a cycle and hence irrelevant. 

Later, Rosenthal developed a rigorous program to regularize scalar and gravitational |12l [53] perturbations 
through second order in perturbation theory. Using these regularized perturbations, he derived formal expressions for 
the second-order gravitational self force in Ref. [25^ in a gauge different from the standard Lorenz gauge, which he called 
the Fermi gauge |23]- The Fermi gauge is a coordinate system in which the regular part of the field perturbation and 
its covariant derivative are made to vanish on the worldline. As such, the first-order self-force corrections vanish in this 
gauge and the background spacetime is composed of the original background plus the first-order metric perturbations 
generated by the motion of the SCO. Unfortunately, Rosenthal's approach does not seem practical for numerical self 
force computations for EMRIs since constructing the Fermi gauge requires the first-order self force on the SCO to 
be made to vanish, which is accomplished by integrating the MiSaTaQuWa equation - a feat that has yet to be 
accomplished. 

Recently^, Detweiler j26J has shown with matched asymptotic expansions that the motion of the SCO can be 
described by a geodesic in a perturbed background spacetime having a smooth metric at the location of the particle 
through second order in e. Detweiler's results imply that one may continue to interpret, at higher orders in e, the 
SCO's motion as either being perturbed by self force corrections on a fixed background or as being geodesic in a 
perturbed spacetime. 



F. Organization 

This paper is organized as follows. In Sectionjlljwe give a brief overview of the effective field theory (EFT) approach 
as it applies to EMRIs. The EFT framework provides an efficient way to systematically calculate the self force at 
higher orders in e. In Section |III| we discuss how to consistently implement outgoing boundary conditions for the 
field within a variational principle. In Section [IV| we demonstrate how to isolate and evaluate the singular part of the 
well-known linear scalar self force on a charge in a non-vacuum spacetime so that we may compare the result from 
our formalism and methods to the standard result in Ref. |27| . In Section |v] we develop a class of nonlinear scalar 
models that is designed to have a structure analogous to the perturbation theory used to describe EMRIs in General 
Relativity. In Section |VI| we then compute the formally divergent self force expressions in the nonlinear scalar model 
through third order in e. In Section [VII| we renormalize these expressions by introducing counter terms into the action 
to cancel those divergences. In Section [Villi we write down the finite, third-order self force equations of motion. In 



Section IX we conclude with a discussion. The Appendices are devoted to deriving the quasilocal expansions used in 
regularizing the self force expressions in Section |VH to proving that power-divergent integrals vanish in dimensional 
regularization and to listing the Feynman rules for the nonlinear scalar model introduced in Section W\ 



The regular part of the self force in the nonlinear scalar model valid through 0{e ) is given in (7.7| and (8.3 1 




using the Detweiler- Whiting decomposition [H] for the retarded Green's function and in ( |8.5| and ( |8.9[ ) using the 
Hadamard decomposition |2H]. Using the Detweiler- Whiting decomposition we find evidence suggesting that the self 
force through 0{£^) can be written solely in terms of the regular part of the field and its derivatives when evaluated 
on the worldline. In a later paper in this series we will explicitly show that this is indeed the case by calculating the 
radiative scalar perturbations and computing their effect on the SCO's motion [29'. 

We use units where c = 1 and define the gravitational constant G in terms of a mass parameter (mpi) to be 
327rG = rripi^. The metric signature is (— , +, +, +). We frequently use the notation where the worldline coordinates 

at proper times r and r' are denoted by the shorthand and z'^ , respectively, so that z^^ = z'^{t) and z'^ = z^{t'). 
The same goes for tensors evaluated on the worldline at some proper time. We also use a mixed notation where a 
quantity such as V{x] z] indicates that F is a function of x'^ but is a functional of z'^{t). 



II. EFFECTIVE FIELD THEORY APPROACH 



In this paper we use the effective field theory (EFT) approach [30] to calculate the self force equations of motion 
for a SCO interacting with a scalar field in a background curved spacetime. The SCO is described by a point particle 
on a worldline with coordinates 2:^(r). 

As a description of gravitational EMRIs, one may worry about the applicability of the point particle approximation, 
particularly its validity for higher-order self force calculations, since it is known that point particle solutions do not 
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exist in General Relativity fSP. The EFT approach is consistent with this statement since the solutions for the metric 
perturbations generated by a point particle approximation of the SCO's motion are not applicable everywhere in 
the spacetime but are valid when the length scales being probed do not reach or smaller. In those cases, the 
EFT must be matched onto a more complete description for the internal structure of the SCO. However, for length 
scales larger than i?„i the point particle approximation is perfectly valid. In fact, the finite extent of the SCO can be 
parameterized in the EFT by including in the point particle action all terms that are consistent with the underlying 
symmetries: general coordinate invariance, reparametrization invariance of the worldline and local 5*0(3) rotations (if 
the SCO is otherwise spherical when removed from all external influences). In this way, our ignorance of the SCO's 
internal structure is merely parameterized and the values of these new parameters are determined through matching 
calculations [30]. In other words, the EFT description of the SCO dynamics essentially utilizes matched asymptotic 
expansions at the level of the action instead of the equations of motion, which is discussed further in Ref. |32] . 

Using EFT methods, it was shown in Ref. [33J that moments induced on the SCO as it moves through the background 
spacetime exert a (tidal) force beginning at fourth order in e if the SCO is a black hole or neutron star. For a white 
dwarf, these finite size effects can be enhanced since a white dwarf is significantly larger than its gravitational radius. 
Hence, a white dwarf is more susceptible to tidal effects that, in turn, depend on how strongly curved the background 
spacetime is at the SCO's location. It was also shown in Ref. [33 that a white dwarf undergoing some form of tidal 
disruption (i.e., Roche lobe overflow or tidal disintegration) may experience a tidal force that is numerically the same 
order of magnitude as a second-order self force correction. Furthermore, based on simple scaling arguments, the tidal 
disruption was found to occur outside the horizon of a supermassive black hole with mass M when M ^ 1.4 x lO~^Af0 
(assuming the mass of the white dwarf is about IAMq), which agrees within a factor of two with the standard lower 
bound of approximately 3 x 10~^ Mq (see e.g., Refs jMl 155]). 

The EFT approach, as it applies to EMRIs, is introduced in Ref. |33] and extends the original work of Goldberger 
and Rothstein |30| who applied it to the post-Newtonian (PN) approximation for slow moving binary sources of 
gravitational waves. In the PN context, a significant amount of work has been done to calculate the conservative 
equations of motion for non-spinning 1301136] and spinning compact objects |37fl45] . The formalism has been extended 
in Ref. to incorporate radiation reaction for generic orbits and to compute the waveforms that will be measured in 
gravitational wave interferometers; see also [47 , 48J for the radiative multipole moments through 3PN for quasi-circular 
orbits with spinning compact objects. In addition, dissipative effects from the absorption of gravitational waves by 
the compact objects can also be treated in the EFT framework |49]l50]. See also Refs. |51|f52] for applications beyond 
the standard two body problem in General Relativity. 

Outside of the gravitational two-body context, the EFT approach has been used to derive the finite-size corrections 
to the radiation reaction on a charged object in Ref. [53 , to study cosmological perturbations in Ref. [IJ, to classify 
horizon geometries and topologies of higher dimensional black holes in Ref. [55 , and to study the Casimir interactions 
of objects in fluid surfaces and interfaces in Ref. [56 , among other things. 

The practical details of implementing the EFT approach in the context of EMRIs (Feynman diagrams, renormal- 
ization, etc.) will be given throughout the rest of this paper. 

III. A CONSISTENT ACTION PRINCIPLE FOR OPEN CLASSICAL SYSTEMS 

The EFT framework utilizes a perturbative treatment at the level of the action. Traditional approaches to the self 
force problem |5J [HI [ST] [57] operate at the level of the equations of motion. We approach the self force calculation from 
an action principle for two important reasons. First, the subtraction of divergences is most easily handled by including 
counter terms in the action, which is particularly useful at higher orders in perturbation theory as we shall see in 
Section |VII[ Second, an action principle unambiguously determines the appropriate worldline for the SCO. Simply 
put, it is the one that extremizes the action. In particular, at first order in e the worldline describing the motion of 
the compact object derived using the EFT approach [53] is given by the self-consistent solution to the MiSaTaQuWa 
equation [S] [S]. However, one may prefer to additionally expand the worldline in s so that the first-order self force 
corrections nudge the SCO away from the leading order geodesic motion. Such an approach is only valid until the 
so-called dephasing time when the corrections induce sufficiently large deviations from the original geodesic that the 
perturbation theory then breaks down. This has caused some confusion in the community about which worldline 
the self force corrects but it is clear using an action principle that the appropriate worldline is the one determined 
self-consistently with the self force, not the perturbed geodesic, which is also advocated in Ref. |58| . 

In this section, we show how to consistently implement retarded boundary conditions when integrating out the field 
from the action. To provide a context, we derive the formally divergent self force in a linear scalar theory using the 
EFT approach. However, to make some of the details and issues more transparent we will not use Feynman diagrams 
in this section but instead calculate the self force step by step. We first attempt to derive the self force from the usual 
action principle, which ultimately fails to describe dissipation in the particle's dynamics, viz., the inspiral evolution. 
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We then construct a new action principle that is capable of consistently and correctly incorporating the field's outgoing 
boundary conditions with the SCO's dynamics. This formalism has been discussed previously in |331 146| within the 
broader language of quantum field theory. Here, we motivate and develop this new action principle using only classical 
arguments. 



A. The usual action principle 



Consider a linear scalar field coupled linearly to a scalar charged particle with mass m and charge g in a non- 
vacuum background spacetime {R^n 7^ 0) so that the action for this system is given by 

Siz^'A] = -^^(0,a0'"+ei?02)-m j dT + q jv{x;z\^ (3.1) 

where = J d^x g^l^ is the integral with the invariant volume element. ^ is a coupling constant and 



V(x- z]= I dT ^ ^ " . (3.2) 



To calculate the self force we first find the wave equation for the field by extremizing (3.1 1 with respect to (j) only, 

a4,~^R(l}^ -qV{x;z]. (3.3) 
The solution with outgoing (retarded) boundary conditions is^ 

<j){x) = q f D,,tix,x')V{x';z] (3.4) 



where D^ctix, x') is the retarded Green's function (or retarded propagator), which satisfies the following inhomogeneous 
wave equation, 

nD,,t{x,x')-^R{x)D,,t{x,x') = ^--^^ ^ . (3.5) 



Substituting the solution in (3.4| back into the action (3.11, which is what is meant by "integrating out" the field in 
this context, gives the effective action |30l I32| . 

S,s[z^']^-m J d-^+Y fj^ Vix;z]D,,tix,x')Vix';z] . (3.6) 

Incidentally, it is this step (integrating out the field) that is performed indirectly by Feynman diagrams. Drawing 
the relevant Feynman diagrams at a given order in perturbation theory amounts to perturbatively solving the wave 
equation and substituting that solution back into the action. The advantage of using the diagrams is that one does not 
have to explicitly solve the wave equation order by order, which is particularly useful at higher orders in perturbation 
theory. 

Calculating the equations of motion for the worldline (equivalently, the self force) from the effective action follows 
from the usual variational principle, 

and yields, upon recalling the identity D,-ct{x' , x) = £'adv(a^, a;'), 

ma'' = ^ / dT'(a^ + P'"'V,)(Z?ret(2^z^') + ^adv(^^,^^')) (3.8) 



^ We ignore homogeneous solutions throughout. 
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where P'^'^ = g^'^ + u^u'^ and we have used the fact that 



Sz^^ir) +^ gi/2 



(3.9) 



Despite imposing retarded boundary conditions on the solution to the wave equation, the resulting effective action 
yields equations of motion for conservative particle dynamics since the Green's function appearing in (3.8|), Dret+Dadv, 



is time-symmetric. This is simply a consequence of the fact that the full Lagrangian for the field-worldline system is 
time-reversal invariant and, as such, yields a variational principle for t he worldline dynamics that is time-symmetric. 
If we couple the field to an auxiliary source J{x) via a term in (3.11 of the form J<p, then the field generated by 



the motion of the particle can be calculated from the effective action using 



5J{x) 



(3.10) 



which here gives 



4>{x) 



{D,^t{x,x') + D^A^{x,x'))V(x' ; z] 



(3.11) 



Therefore, the effective action yields a radiative field that no longer respects the outgoing boundary conditions and 
disagrees with the explicit solution to the wave equation in (3.4|. In fact, (3.111 describes radiation with no net fiux 



of energy leaving the system, which is consistent with the absence of a dissipative component to the self force on the 
particle in (3.8 1. 



It is worth remarking that for binaries with comparable masses in the post-Newtonian (PN) approximation the 
effective action is just the so-called Fokker action [591 [60] . The Fokker Lagrangian has long been known to describe 
only conservative dynamics (see e.g., Ref. [SI])- In Section IIIB we introduce a new action principle that essentially 
generalizes the Fokker action to allow for a description of dissipative processes as well as conservative ones. We have 
applied this new action principle in Ref. |1H] to derive the 2.5PN radiation reaction force of Burke and Thorne |62 [ [55] 
from a Lagrangian formulation. 



B. A consistent action principle for open systems 



The shortcomings found in the usual action principle, as it applies to open systems, are caused by the time-reversal 
invariance of the full Lagrangian. One way to incorporate time-asymmetry into the true evolution of the system is to 
compare the actions of two different histories of the worldline and field variables, call them (zf,(/>i) and (2:3, (/'2)- A 
useful comparison is to compute the difference in the actions of the two histories 

5[z^4,</)l,(/.2] = 5[zr,0l]-5[z2^<^2] . (3.12) 



The interpretation of (3.121 is as follows. Since the Lagrangian is time-reversal invariant it follows that the second 
term on the right side of (3.12) is 



'2 1 V2\ = 



dtL[z!^,z!^,(f)2,do(l)2 



(3.13) 



Hence, the second history evolves backward through time from its final configuration to its initial one. This is why we 
have formed the difference of the actions - (3.121 naturally compares the forward time-evolution of the first history 



with the backward evolution of the second history. Thus, (3.121 provides a measure of the time-asymmetry between 
the two histories. 



Obviously, if the two histories are the same, then (3.121 vanishes identically - if we compare the history evolved 
forward in time with the same history evolved backward in time then there is no difference. In addition, (3.121 is a 



legitimate action in its own right since its extremal points correspond to extremal points of the actions for the two 
histories. However, there is only one relevant and physical history (for a given set of initial data), which we denote 
by (z'^, (j)). Since the extremal points of the action determine the equations of motion then it is reasonable to require 
that the histories are equal to the physical one after performing the variations so that 



02 = 



(3.14) 
(3.15) 
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In this way, any time-asymmetries between the two histories will be captured in the description of the physical history's 
dynamics. 

Incidentally, the action in (3.121 is the one resulting from the classical limit of the corresponding quantum theory 
(of a quantum scalar field interacting with a point particle source) as derived using the so-called "in-in" formalism. 
The in-in formalism, first introduced by Schwinger [64J and Keldysh [65j, has been proven to give a causal description 
for open quantum systems. Hence, its classical limit provides a consistent implementation of the outgoing radiation 
boundary conditions. See Refs. |33l fiSl [55^ and references therein for more details. 

Let us now follow the same steps as in the previous section to calculate the self force but starting instead from the 
action in (|3.12|. Written out fully (|3.12| is 



5[zf 2, 01,2] ^ - mj dTi+mJ dT2-^ ^(01.a</'l'" + + I ^(02,a02'" + ^R4'l 



+ q V{x;zi] 



V{x;z2]4'2 



(3.16) 



where ti^2 are the proper times associated with the worldlines z^j- It is convenient to make a change of variable to 



6_ = 01 - 02 

^+ = ^(01 + <P2) 



(3.17) 
(3.18) 



in which case the action (3.161 now reads 



m I dri + m dr'^ 



.+,^0_''^ +^i?0+0_)-Kg / V+{x;zi,2]^.+q I Vl(x;zi,2]0+ (3.19) 



where 



The wave equations for (j)± are 



V-{x; zi^2] = V{x; zi] ~ V{x; Z2] 
V+{x;zi^2] = ^{V{x-,zi] + V{x;z2]) ■ 



□ 0± - ^i?0± = -qV±{x; zi,2] 



(3.20) 
(3.21) 

(3.22) 



When reducing the histories to the physical one using (3.141 and ( |3.15 1 it follows that the field in the ± variables 
becomes — > while 0_ — > 0. This implies that (i.e., the average of the two field histories) reduces to the 
physical field in this limit and the solution to the corresponding wave equation above must satisfy retarded boundary 
conditions. Therefore, 



(3.23) 



(j>+{x)^q D,ct{x,x')V+{x';zi^2]- 

J x' 

Solving the 0_ wave equation for momentarily unspecified boundary conditions gives 

0_(2;) = q [ D{x,x')V-{x'-,zi^2] 



(3.24) 



for some Green's function D{x, x') that will be determined shortly. Putting these solutions into the action (3.19 1 gives 
the effective action, 



PL {x; zi,2] {D,^t{x, x') + D{x\ x))V+{x' ; zi ■ 



(3.25) 



The equations of motion for the particle follow by extremizing the effective action with respect to either z^ or Z2 and 
setting z^2 = using (3.151, which gives 



ma'' = |- / dr'(a'' + P'"'V,)(A■et(^^^''')+^(^''^^'')) • 



(3.26) 
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In order for the equations of motion to exhibit the proper causal structure D{x, x') must equal the advanced propagator, 

D{x,x')^D,d,{x,x') (3.27) 
upon using I?adv(2;', x) = Drot(2;, x'). Therefore, the (formally divergent) self forced motion of the particle is given by 



/ dr'(a^ + P'^''Vji?,et(^^,^^') 



(3.28) 



We emphasize that what we have done here is to integrate out the field in the new action (3.12 1 to obtain an effective 
action whose extremum describes particle motion with dissipative effects. 

The fact that (/)_ satisfies advanced boundary conditions (according to (3.24 1 ) is not a problem because the source of 
the field vanishes when setting — zi^ — implying that itself vanishes. Therefore, there is no contribution 
to any dynamics from advanced boundary conditions. For completeness, the effective action is 



'S'cff[^:f,2] = -™ JdTi+m JdT2 + q'^ J J V^{x; Zi^2]Drct{x,x')V+{x'; Zi^2] 



(3.29) 



Note the similarities and differences between (3.29) and (3.6 1. 

Incidentally, coupling the fields and (j)2 to external currents Ji and J2, respectively, computing either of (55eff/(5Ji,2 
and then setting = Ji,2 — gives the field radiated by the particle 



0(^) 



Drctix,x')V{x';z] , 



(3.30) 



which agrees with (3.23) when applying (3.14) and (3.15) and demonstrates the internal consistency of the formalism. 



C. Formal developments 



Here, we briefly introduce some notation for the formalism developed in the previous section that will be used in 
the following sections to compute higher-order self force corrections in a nonlinear scalar model of EMRIs. 

The action in (3.16) can be written in a condensed form by introducing a "metric" cab used to raise and lower the 

(3.31) 



1, 2 indices labeling the histories. 



CAB 



1 

-1 



„AB 



where a capital Roman letter takes values in {1,2}. Then, (3.16) can be written as 

'5'[<2>i,2] ^ -m J dTi+m J dT2-^ J i(l}^^(t>'X + ^R<P'^(t)A) + qJ^V'^{x; z\(j)A 



(3.32) 



with Va{x\ z] = V{x; za], (j)^ = 01, = —4>2 and repeated indices indicate summation as usual so that = c^^Vb 



and (j) 



nAB, 



>B- Likewise, using the "metric" 



Cab 



1 

1 



(3.33) 



where a lowercase Roman letter takes values in {+, — }, allows for (3.19) to be expressed more compactly as 



'5[<2> ^1,2] = -m J dn+mj dT2-^ j}"^"^^"^'^ ^R<l^"4>a) + qjv\x- z]<j>a (3.34) 



with V-=Vi- V2, V+ = [Vi + ^2)72, (/)+ 



and 6 — 



the indices that label the varia 
observing that 



Notice that (3.32) and (3.34) have the same form, suggesting that the action (3.12) is, in a sense, "covariant" in 

ales. The transformation A that goes from the ± basis to the {1,2} basis is found by 

(3.35) 



^ 2 
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or (pA = ^A"'4'a where can be read off from the equation above. In addition, the metrics cab and Cab are related 
by A through 

0,0,0"" = 20+0_ =cl>l-cf,l = Mbc^"" = (AA"0a)(ABV6)c'^'' (3.36) 

from which it follows that 

as could have been guessed from the "covariant" structure of the action with respect to the histories' labels. 

The ± basis gives a convenient representation for the effective action in (3.291. To see this, define the following 
matrix of propagators 



Di,dvix,x') 

Drctix,x') 



(3.38) 



which is symmetric under interchanges of indices and variables so t hat D °'^{x, x') = D^°'{x\ x) upon using the identity 
Da.dfW ^x) = D^ctix,x'). [Here, D ^{x,x') — D,-ct{x, x').\ Then, (3.25) can be written in a more symmetrical form 



Scfi [^1 ,] 



Va{x;z]D''''{x,x')Vb{x';z] 



(3.39) 



IV. SELF FORCE REGULARIZATION IN A LINEAR SCALAR THEORY 



The self force expression in (3.28 1 is formally divergent when the retarded propagator is evaluated at t' = r. Hence, 



a suitable regularization procedure must be given to ensure a sensible and regular self force on the particle. In this 
section, we describe how we regularize the self force. We use the Hadamard decomposition of the retarded propagator 
into its so-called direct and tail pieces |i28_ and then proceed to evaluate the singular parts of the proper time integral in 
(3.281 using quasilocal expansions for r' near r (see Appendix [A| . We demonstrate that our regularization procedure 
yields the correct finite self force expression as first derived by Quinn in [57] . When calculating higher order self force 
corrections in Section VI we use the Detweiler- Whiting decomposition [T3] of the retarded Green's function as this 
has better behaved properties than Hadamard's form. The impatient reader who wants to get on with the nonlinear 
scalar model may safely skip ahead to Section \V\ 

Whenever x'^ and x''^ can be connected by a unique geodesic (the set of all such points connected to constitutes 
the normal neighborhood of x^^), the retarded propagator can be expressed as |281 f57| 



D,,t{x, x') = j^O+{x, J:,,)(^A^/\x, x')d{a{x, x')) + V{x, x')B{-a{x, x'))^ , 



(4.1) 



which is the sum of a "direct" part (first term) and a "tail" part (second term). Here, Yi^' represents a space-like 
hypersurface containing the point x' , 9{x, S^;/) equals one if x'^ is to the future of E^;/ and zero otherwise, A{x, x') is 
the van Vleck determinant, a{x, x') is half the spacetime interval along the geodesic connecting the two points (also 
called Synge's world function) and V{x^x') is a regular function. For notational convenience, let 



DA,r{x,x') = ^e+{x,i:.,,)A^'^{x,x')5{a{x,x')) 
Ait 



(4.2) 



be the direct part of the retarded propagator, which only has support on the future null cone of x'^ . 

If Tin and Tout represent the proper times at which the worldline enters and exits, respectively, the normal neigh- 
borhood of the point a;^ = z''(t) then we can split the integral in (3.281 into contributions from outside and inside 
the normal neighborhood of z^{t), 



dr' 



dr' 



dr' 



(4.3) 



For the integration inside the normal neighborhood we may substitute Dj-ot from (4.1) so that (4.3) becomes 

(•Tout 



DiUz'',zn + Sir-T')V{z^',z'^) 



(4.4) 
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where we have used 9{—a{z^, ) = 1 for all r'. The contribution proportional to V is regular and can be combined 
with the second line in (4.4 1, which is also regular, so that 



c— >0+ , 



(4.5) 



where the limit is used to ensure that the second line receives no contribution from the singular coincidence limit of 



the propagator. All of the divergent structure is contained in the first integral of (4.5 1, which we evaluate next 



There are two potentially divergent integrals that appear in (4.5), which are proportional to 



/oo poo 
dr' DdUz'',z''') and P'"' dr' V ^D^^^iz" , z"' ) 
~00 J —OO 



(4.6) 



Each integral has support only when t' = t (which is why we have replaced the integration interval from [Ti„,Tout] 
to the whole real line) so that we can expand the integrands about s = t' — r = 0. Using the expansions listed in 
Appendix \K\ it follows that 



/oo 1 /"^ 

dr' Dd,,{z'',z''') = — ds 
-oo 47r 



Sis) 



A 

4^' 



(4.7) 



which is completely divergent since there are no 0(5") terms in the product of A^/^5(cr). 

The second divergent integral in (4.6) can be calculated in a similar way. Using the relations in Appendix [A| we 
find that 



dssgn{s)6is) . (4.8) 



The integrand in the second term is proportional to a product of distributions, sgn(s) and S{s), suggesting that 
the integral is not well-defined. Specifically, the delta function has support only at s = 0, which is where the sgn 
distribution is not defined [55] ■ We will return to this issue shortly. 



The regularized self force follows by substituting (|4.7| and (|4.8|) into (|4.5| giving 

^2 



fa^' ( A 



1 



I Da, 



A^J-i^'-'^U^-^^'^" ' 3dr 



g2 lim / dr'(a^ + P^''V,)Act(z^,z''') 







ds sgn(s) d{s) 



(4.9) 



To evaluate the remaining s-integrals, we remark that the divergence arises when s = or, equivalently, t' = r. It 
is natural to parameterize the limit r' — > r by shifting the argument of the delta function in (4.9 1 from s to s + T 



for some T ^ 0+. This prescription also resolves the problem associated with (4.8) since now the delta function has 



support at s = — T where the sgn distribution is unambiguous. The two integrals in (|4.9| then evaluate to 

S{s) 



A = 



ds 



lim 



" ds'±±n. Hm 1 

S T^0+ T 



/O 
ds sgn(s) (5(s) — > Hm / ds sgn(s) 5(s + T) = — 1 
-00 T-^o+ 



(4.10) 
(4.11) 



and the self force becomes 
1 



gMim / dr'(a^ + P^'^V,) Act z^') . (4.12) 



The divergence in the self force is proportional to the acceleration. Hence, by introducing mass counter terms into 
the effective action (3.19) of the form 



(4.13) 
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we find that their contribution to the self force is —SmO.^ and we can render (4.121 finite if 

" 2 47r " 8^ T "o+ T ' 



(4.14) 



If instead of using the prescription s — ^ s + T we used dimensional regularization, we would find that the divergent 
integral in (4.7) actually vanishes so S„i is zero. As is well known in the field of high-energy physics, an integral that 
diverges as a power of a cutoff (here, T is the cutoff parameter) will vanish in dimensional regularization^. For the 
sake of completeness, we discuss in Appendix [B] how to use dimensional regularization in this context and verify our 



claim that ( |4.7[ ) vanishes. 

Putting the pieces all together yields the (finite) self force equations of motion for the particle, 



6 47r 3 47r dr 



qMim / dT'(a'^ + P^''V,)A■et(^^2^'), (4.15) 



which is the expression first derived by Quinn in [27,. From the equations of motion, one can define an effective mass 
for the particle as 



Weff (t) — m — q lim 



dr' Aot(2^^'^') 



from which the rate of change in the particle's mass is 
drucs 



dr 



A|:(l-60i?(^^)~ A" lim 
12 47r e->-o+ 



dr' V^D,,,{z^',z''') 



(4.16) 



(4.17) 



See Ref. [37] for a general discussion and Ref. [69J for an example of the scalar charge's mass evolving in a cosmological 
context. 

In the sections below, we will determine the higher order corrections to (4.15)-(4.17l but we must first generalize 
the action in (3.1 1 and (3.12) to describe a nonlinear field theory model appropriate for extreme mass ratio inspirals. 



V. A NONLINEAR SCALAR MODEL OF EMRIS 



There are many nonlinear scalar field theories that admit higher-order self force effects. However, we are interested in 
a class of models that possesses a structure analogous to the General Relativistic description of EMRIs. In particular, 
the motion of a SCO about a supermassive black hole is described by the perturbed Einstein-Hilbert Lagrangian (in 
the Lorenz gauge for the trace-reversed perturbations) about a black hole background spacetime (with metric g^i, and 
dimensionless perturbations h^^/nipi) and is given schematically by^ 



(5.1) 



where m^j = l/(327rG) is a mass scale associated with the gravitational constant G and the coefficients a"'^^^ "^^" (x) 
are dimensionless and depend only on the vacuum background metric. 

One class of nonlinear scalar models that admits a structure analogous to (5.1) is 

oo „ 

= - E / af{^)dacl^dpcj,r-' (5.2) 



where a"^ = ang"^ is proportional to the metric since that is the only dimensionless rank-2 tensor available to contract 
with the derivatives on the fields and we take a2 — 1- 



^ The field that generates the divergence must be massless also. 
We have explicitly verified the structure of these higher order terms in the Einstein-Hilbert action through sixth order in h/^i,, which is 
more than sufficient for our third order calculations here. 
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The action of a point particle coupled to gravitational perturbations in General Relativity is given by expanding 
the integral of the worldline's proper time on the perturbed spacetime, 



(5.3) 



where the (r) depend only on the particle's four-velocity u"(t) and 6q = 1. In the EFT approach, extra terms 

describing the effects from the finite-size of the SCO can be included [30 . However, we will leave the discussion of 
finite-size effects for a future paper. The analogous action for the interaction of the particle with the nonlinear scalar 
field is 



oo „ 
n=0 ' Pl 



(5.4) 



where the 6„ are constants and we take 6o = 1. 



Combining (|5.2| and (|5.4| gives the full action for an analogous nonlinear scalar model of gravitational EMRIs 

0,Q0'"A^((?!)/mp;) - m / dr B {(t) / nipi) 



where 



2a. 



n+2 



^1 + \mpi 



B{cf>/mpi) = l + J2 



The wave equation for cj) and the corresponding equations of motion for the particle are 



U(h = -- 



A' 



A mpi rupi 
a'' = -P^'-'V^lnS 



dr 



S^ixt" - z^(t)) B' 



jl/2 



(5.5) 

(5.6) 
(5.7) 



(5.8) 
(5.9) 



where a prime denotes differentiation with respect to the function's argument, viz., tp/nipi. 



Field redefinition 



The action in (5.51 can be greatly simplified through the following field redefinition. Let '4>{x) be given by 

V'.Q = 4),aA{(l)/mpi) (5.10) 

or -0 = F{(t>/mpi) + C where F(x) = J dx A{x) is the anti-derivative of A and C is an int egration constant that fixes 
^ = when cj) — 0. If A cannot be integrated in closed form then one can integrate (5.10 1 to find -0 perturbatively in 
terms of </> using (5.6) 



rUpi rUpi 12 \mpi 



3a4 — 

216 \mpi 



and vice versa, 



mpi iripi 12 \mpi J 216 \mpi 



The change of variables implies that the action in (5.2) becomes 

1 



(5.11) 



(5.12) 



(5.13) 
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The key point to note is that while we constructed (5.2 1 to mimic the nonHnear structure of perturbed General 
Relativity we have found that (5.2 1 can be transformed to a linear field theory, that is, one where the field does not 
interact with itself. This will have important implications regarding the divergent structure of higher-order self force 
effects in this model (see Section VII I and the corresponding renormalization of coupling constants. 

A curious question naturally arises: If the nonlinear scalar theory is analogous to perturbed General Relativity then 
to what extent can one make a field redefinition to remove some, if not all, of the higher order (n > 2) terms in (5.1 1? 
We will address this in a forthcoming paper. 

The full nonlinear scalar model for EMRIs is described by the following action upon using the change of variables 
from <f> io tp, 



where 



n! \mpij 



and the first three {c„} coefficients in terms of the original coefficients, {a„} and {bn }, are 

ci = bi 



1 



asbi 



1 2 

^3^2 - TTia^bi + -a^bi . 



C2 = b2 
C3 = ^3 

The corresponding equations of motion for the field -ip and the worldline are 

(54(j;^ - 0^(t)) 



1 

12' 



rripi 



dr- 
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1/2 



-C'iiP/mpi) 



where, as before, a prime denotes differentiation with respect to the argument of the function, viz., ip/nipi. 



(5.14) 
(5.15) 

(5.16) 
(5.17) 

(5.18) 

(5.19) 
(5.20) 



B. Power counting and Feynman diagrams 



The relevant scales in this nonlinear scalar model of EMRIs are the size of the SCO Rm, the curvature scale of the 
background spacetime TZ and the gravitational constant m^;. These scales form a natural dimensionless parameter 
e = Rm/TZ that can be used to form a perturbative description for the EMRI whenever <C TZ- If the SCO, 
with mass m, is in the strong-field region of a supermassive black hole, with mass M, then e is just the mass ratio 
m/M. Calculating the self force in the EFT approach using Feynman diagrams requires knowing which diagrams will 
appear at a given order in e. This is accomplished by power counting, which determines the dimensional scaling of 
the interactions in (5.141 in terms of the relevant scales in the theory. 

The field ip{x) describes the radiation from the SCO implying that daip/ip ^ l/TZ and hence ^ TZ. Requiring 
that the kinetic term for the field is the leading order term in the perturbation theory (i.e., /^(Si/))^ ^ 0(e°)) implies 
that i/j ^ l/TZ. The worldline integrals in (5.141 then scale as 



1 



(5.21) 



For n — this is —m J dr ^ mTZ =: L and describes the leading order scaling for the SCO dynamics since the 
interactions with the field are perturbative corrections to this. In terms of e = Rm/TZ ^ ■m/{rn^{TZ), (5.21 1 becomes 



mcn 



dri^'^z") - L 



n/2 



(5.22) 



In Appendix [C] we give the Feynman rules for the theory described by the action in (5.141. The vertices appearing 
in Feynman diagrams encode the interactions between the field and the SCO and are described by the worldline 
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(a) Loop diagram 



(b) Disconnected diagram 



FIG. 1: A loop diagram (a) and a disconnected diagram (b), neither of which can be generated by a local classical theory. The 
solid line represents the worldline of the small compact object and a curly line represents the field's propagator. 



integrals in (5.141. which scale with e as shown in (5.221. Therefore, to compute the effective action (and hence the 



self force) at 0(e^) requires assembling diagrams with certain combinations of vertices such that the entire diagram 
scales as . These vertices are then connected by curly lines denoting the appropriate Green's function for the field 



(recall that there are two appearing in the matrix of propagators in (3.38)). An example is shown in Figure [2] for the 
Feynman diagram that contributes at first order in e to the effective action. Notice that, in this figure, the labels a, b 
refer to the histories, discussed in Section [ill B[ 

There are many possible diagrams one may draw. Since Feynman diagrams correspond to terms in the effective 
action, which are computed by perturbatively solving the wave equation and substituting that solution back into the 
full action, then the diagrams must be drawn according to that procedure. The solutions to the wave equation will 
never generate terms in the effective action involving ~ J dr D{z^ , z^), for example. Such terms would appear as 
loops of curly lines (see Figure l(a)[ ) and are not admissible in a classical theory such as (5.14). In fact, diagrams 
with such loops of curly lines turn out to be quantum effects and are suppressed in the effective action by powers of 
h/L, as discussed in Therefore, all diagrams here must not contain any loops of curly lines. 

In addition, all diagrams must be connected in the sense that cutting the worldline at any place leaves the resulting 
diagram connected by curly lines. Classical solutions to the wave equation cannot generate terms in the effective action 
corresponding to disconnected diagrams. To see this, let's assume they did. Then there could be terms appearing in 
the effective action like the one in Figure |l(b)| and schematically written out here as 

-I 2 



dr / dr' Diz^'^z^' ) 



that could only come from the following term in the full action 

n 2 



dT^ljiz^") 



^ dr dr' iP{z^')iP{z^' ) 



(5.23) 



(5.24) 



Clearly, such terms are not present in (5.14) since the Lagrangian is local in time. Therefore, such terms also cannot 



appear in the effective action and the relevant Feynman diagrams must be connected ones. 

These two conditions together (no loops or disconnected diagrams) imply that every diagram contributing to the 
effective action scales as a single power oi L = mTZ. Essentially, these rules come about because (5.14) describes a 
local, classical field theory. 



VI. SELF FORCE IN THE NONLINEAR SCALAR MODEL 



In this section we calculate the expressions for the self force on the compact object through third order in e within 
the nonlinear scalar model of (5.14), i.e., after implementing the field redefinition. The Feynman diagrams that arise 



at each order in this perturbation theory have the same topology as a subset of those diagrams appearing in the 
gravitational case describing the dynamics of physical EMRIs. As a result, the computational tools and methods used 
and developed here will be indispensable for calculating the gravitational self force at higher orders in the Lorenz 
gauge in a future paper. 

To show that our renormalization procedure works for singularities arising from nonlinear interactions of point 
sources we will explicitly retain all divergences until we are finished computing the third-order self force. At that 
point we will shall identify the counter terms that cancel these divergences and render the self force expressions finite. 



Self force at first order 



At first order there is a single Feynman diagram contributing to the effective action, which is shown in Figure 
[2] We have already regularized the first-order self force in Section |IV[ In that section we decomposed the retarded 
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Cl 



Cl 



FIG. 2: The Feynman diagram contributing to the effective action at first order in e. Power counting indicates that this 
diagram scales as ~ eL. 



propagator into its direct and tail pieces in order to obtain a sufficiently regular self force expression to compare with 
the known result in |27| . However, that decomposition is not particularly useful for practical self force computations, 
especially at higher orders, as we shall see in Section [VlII| 

A better decomposition is provided by the Detweiler- Whiting prescription [T^ where the retarded propagator is 
split into a regular (i?) and singular (S) part. 



Drct{x,x') = Dr{x,x') + Ds{x,x') 



(6.1) 



The advantage of this decomposition is that the regular part Dji satisfies the homogeneous wave equation, is regular 
on the worldline and is the part of Z^i-ct actually responsible for forcing the SCO. Therefore, the corresponding 
contribution to the self force is regular on the worldline and one can make sensible predictions. The singular part Ds 
satisfies the inhomogeneous wave equation, carries all the divergent structure of D^^t ^^nd exerts no force on the SCO. 
When X and x' are in a normal neighborhood we can write the regular and singular parts as 



(x, x') = ^Ai/2(x, x')5{a) [e+{x, - e^{x, E,0] + jz'^i^^ 



e+ix,j:.,:)ei-a) + ^eia) 



{x,x') = ^A^^^{x,x')S{a) - ^V{x,x')0{a) 



Stt 



Stt 



(6.2) 
(6.3) 



See [S3 for a pedagogical presentation and discussion of the Detweiler- Whiting decomposition. 

Following steps similar to those implemented in Section IV (but using the Detweiler- Whiting decomposition instead) , 
we isolate the singular and regular parts of the integrals in (3.281 so that 



/OO 
dT'Dn{x'',zy) 
-oo 

and, by construction, is regular on the worldline so that V i,Ir{x'^) is well-defined when x^ 
first-order self force is then 



where Ir is defined as 



a" -t- 



2 ? 



"pi 



(6.4) 
(6.5) 

(6.6) 

[12] . The regularized 
(6.7) 



Notice that (6.7 1 agrees with our earlier derivation in Section IV based on Hadamard's decomposition of the retarded 



Green's function, which is expressed in terms of direct and tail pieces. To see this, we identify q = mci/mpi so that 
the tail part of the self force in (4.121 is related to the Detweiler- Whiting regular part by 



(a^-|-P^'^V,)/fl(z^ 



lim 

£-^0 + 



dr' (a'^ +P^''V,)Aet(^'^,z'^') 



B. Self force at second order 



At second order in e, power counting indicates that one diagram contributes to the effective action, which is shown 
in Figure [3] and scales as e^L. The Feynman rules for that diagram imply that the corresponding terms in the effective 
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FIG. 3: The diagram contributing to the eflective action at second order in e. Power counting indicates that this diagram 
scales as ~ e^L. 



action are 



2! 



X J x' J X 



Taix; z]D''\x, x')Tbc{x'; z\D'%x' , x")Td{x"; z 



(6.9) 



Here, we use the notation from Appendix[C]so that a vertex with n scalar fields attached to the worldline is represented 

by 



(6.10) 



where Cai-a,," = l^a^"^^ •••Aa„+/"A^BCAi-A„'' and 



1 ^i = ...= 

-1 Ai = ...= 
otherwise 



Ap — 2 



(6.11) 



Substituting in (6.10) into (6.9 1 gives 

2 



mcl 
rripi 



mc2 



Va[x- z]D''\x, x')cbc'V,{x'- z]D'\x\ x")Vd{x"- z 



We perform the contractions in the ± basis and use the non-zero elements of the c-tensors, 



c+-^ = c+— = 1 
C-++ = = 1 

C ^ = C |- = 1 

C++" = C+++ = 1/4 



(6.12) 



(6.13) 
(6.14) 
(6.15) 
(6.16) 



to find that the effective action becomes 



^ 9 

m c^C2 



X 'J X 'J X 



2V^{x-z\ D,,t ix,x')V+{x';z] D,,t {x' ,x")V+{x" ■ z] 
+ V- (x: z]D,,tix, x')V+ix'; z]D,,t{x, x")V+ix": z] + 0{V^) 



(6.17) 



Th e sel f force at second order in e is then found by v arying (6.171 with respect to either or Z2 using the identity 
in (3.9 1 and then setting z^ 2 — z^^ according to (3.151, 



2m„, ^ 



s 2 

TO CJC2 



dr' D,,,{z^,z^') 



(6.18) 



The 0{V^) terms in (6.171 do not contribute to the equations of motion since the variational derivative of those term s 
with respect to or z^Ts proportional to l/_ (x; z] (x z^ — Z2, which vanishes when we set zf = Z2 = z^ from (3.15 1. 
It will be convenient to evaluate each line of (6.18) separately, which we do in the next subsections. 
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1. Regularization of F^^^s^ 



Here, we regularize the divergent integrals appearing in the first line of (6.18 1, 



2m' 



pi 



dT'D,,t{z'',z''') 



The regularized expression for F^^^^ follows simply from substituting (6.4 1 and (6.5 1 into F^^^y which gives 

m?'c\c2 ( A a'' 



47r 2 



(6.19) 



(6.20) 



2. Regularization of F^^^-^ 



Turn now to regularizing the integrals in the second line of (6.181, 



^ 9 

m c^C2 



(6.21) 



Notice that P(26) involves the divergent integral of a divergent integral as opposed to products of divergent integrals 
as in F^^^y One way to deal with these integrals is to insert a factor of 1 = / dfS{f — r') into (6.21 1 to "disentangle" 
the r' and r" integrals. 



^ 9 

rrf^CiC2 



ml 



dfla^i / dT'5(T'-f)A■et(^^-'^') dr" D^.tiz^z'^") 



+ [P^'' dT'<5(T'-f)V,Act(^^^'^') / dr"Aet(2^^'^") 



(6.22) 

The integrals involving the factor of 5{t' — f) can be integrated by letting t' = r + s and expanding the integrands 
about s = 0, as before. Doing this, we find 



dr' 5{t' - f) Act(2^ z^') = ^6{t - f) + Dr{z^, z'^) 

47r 



dr' 5{t' -f)V,D,,t{z^,z^' ) 



A ai, 
4^^ 



<5(T-f) + V,i?fl(z^^^') 



(6.23) 
(6.24) 



Inserting these expressions into P^^^) and integrating over r then gives 



m^CiC2 ( f A \ fl' 

m' 



(il)? / ^^"A■°t(;^^;^''") + K + ^^'^V,) y dT'i?fl(z^z^') j dr" D,,,{z^^\z^^") 



(6.25) 



The remaining integrals over r" have been computed already in Section VIA Substituting in their values from ( |6.4[ ) 
(6.5) and collecting terms yields 



^^c\c2 \ f A Va'' 



47r/ 2 V47r/ 2 



{zn 



+ (a^+P^^V,) y dT'i^fl,(z'',z^')/^(z^')| 



(6.26) 



For the third order calculations below, it is useful to collect an intermediate result that can be inferred from the above 
manipulations, namely 

2 

dr'dr" D,,t{z''.z^^')D,,t{z^\z''") 



A 

in 



P^"- / dr'dr"V,Aet(^^,^^)Aet(^^,;2^ ) 



1 / A 



2V47r 



2(^^yn{zn + J dT'DR{z^,z^')In{z'^') (6.27) 

Rizn 



1 / A 



2 V 47r 



P^'^V.y dT'DR{z^^,z^^')In{z^^') . 



A)p-v./. 



(6.28) 
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FIG. 4: The diagrams contributing to the effective action at third order in e. Power counting indicates that these diagrams 
each scale as ~ e^L. 



C. Self force at third order 



At third order in e there are two diagrams contributing to the effective action and are shown in Figure |4j The 
Feynman rules imply that the corresponding effective action is 



eff 1^1,21 — 21 



Ta{x; z]D''\x, x')TUx')D-''{x\ x")Tae{x")D''f {x'\ x"')Tf{x"') (6.29) 

Ta{x- z]D''^{x, x')ncd{x')D''{x\ x")T,{x")D''f{x', x"')Tf{x"') . (6.30) 
In the ± basis, the nonzero elements of the c^icdg tensor, which appears in the expression for the Tbcd vertex, are 



1 
3! 



X J X' J X" J X 



1 = c + = c__+_ = c_+__ = c+ 

1/4 = C-+++ = C+-++ = c++_+ = c+++_ 
so that the part of the effective action that contributes to the worldline equations of motion is 
o(3)r M 1 _ rn^clcl 



(6.31) 
(6.32) 



X J x' J X" J X 



+ 



^^pl J X J x' J x" J x'" 



V-{x; z] A-ct(a;, x')V+{x'; z\D,,t(x' ,x")V+{x" ■ z\D,,t{x\ x"')V+{x"'- z] 
+ F_(x; 2]D,et(a;, xOl^+(x^ z] Act(a;^ x'Ol^+(a;'^ z]Act(x, x''OV+(a:''^ z] + 0(V^_2) 

3F_ (x; z]L»,et(a;, a;')F+(x') Aet(x', a;'')V^+(a:"; z] Aet(a^', x''')"^^+(x'''; z] 
+ V.{x; z]D,,t{x, x')V+{x'- z]D,,t{x, x")V+{x"; z]D,,t{x, x"')V+ix"'; z] + 0{V^) 



(6.33) 



Variation of (6.331 gives the contribution to the equations of motion from the third order self force, 

^4^2 ^2 



flA'+pM-VJ dr'dT"dT"' D,,tiz'',z'' )Drct{z'' )D,,t{z^ ,z'^ ) 



4 2 2 r r 

+ ""^^(a^ + P^'-V,) j dr'dr'" D,,t{z^,z^^')D,,t{z^,z^"') J dr" i?,et(z^', z^") 



m^CiCs 



a^ + P^-W,) / dT'A-ct(z^,z^) 



dr" A-ct(z'^',z^") 



1 3 



dr' A-ct (z^,z^ ) 



(6.34) 
(6.35) 



We next regularize each line appearing in Fj^^y Only two new integrals need to be regularized and both appear in 



(6.5), (6.271 and (6.28) 



F^^y The remaining integrals in (6.34) can be expressed in terms of the previously regularized integrals in (6.4) 
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1. Regularization of F^^^s^ 

Regularizing F^^^-^ follows similar manipulations and steps as in regularizing ^(26) f'^om Section 
"factorize" the proper time integrals by inserting 1 = / dfS{T' — f), 



VIB2 



in that we 



^(3a) - 



TO, 



J dfi^a^^l^J dr' 5{t' -f)D,,,{z^,z^')^(^j dT"dr"' Aet(;2^ ^^")Act(2^" , ^''"V 

P'"- fdr'Sir' - f)V,A■et(^^^^')) ( /drV Aet(z^^^")Aet(;2^",2^"') ) } . (6.36) 



F, 



II III ^ 



Using (|6.23| and (|6.24|) to regularize the r' integrals gives 

(a^ + F^^'V.) I dr' DRizf^^z^') J dr" dr'" D,,,{z''' , z''" )D,,t{z''" , z^'" ) . (6.37) 



(3a) ~ ^6 



m^:.! I \ 47r / 2 



The remaining integrals over r" and t'" are given by (6.271. Substitution then gives 

2 



"pi 



A 



5 / A 



K^)^''"^''/ ^^'^i?(^"'^'^Vfl(^^')+K+^^''V,) y dT'dT"i^fl,(2^2''>fl(z''',Z^")/^(z''")| . 

(6.38) 



^. Regularization of F^^^^-^ 

Isolating the singular pieces in -F^^f,) turns out to be relatively simple since some of the integrals already factor 
giving 



(36) 



^^l^a^^(^j dr'dr" D,,,{z^,z^')D,,,{z^' ,z^")^(^j dr'" D,,,{z^ , z^^'" )^ 

P^'' j dT'dT"^,D,,t{z^',z^'')D,,t{z^^',z^'")^ (^J dT"'D,,t{z^,z^^"') 
dr'dr" D,,t{z^,z^^')D,,t{z^' ,z^^") \ fp^^ / dr'" V, Aet(^^, ^^"') 



(6.39) 



Substituting the expressions (|6.4|, (|6.5|), (|6.27| and (|6.28|) into F^^^^^^ then gives 



4 2 9 

{3b) ^ 



1/A 

21 47r 



(z^) + 3( i^)p^^''VJnizniRizn 



A 

47r r 



J dr' Z?fl(z^ z^')/fl(z'^') + (^)^^"V, y dr' Dfl(z'', z^VflX^'^') 
+ /,j(2^)(a^ + P^''V,) / dr'i?fl.(z'',z''')/fl(z^') + P^''V,/fl.(z'') / dr' i?fl(z^ z'^')/fl,(z'^' ) \ . (6.40) 
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FIG. 5: Feynman rules for the counter terms that cancel the divergences in the self force through third order in e. Here, S,n is 
a mass counter term while Si and 82 can be associated with renormalizing ci and C2. 



3. Regularization of F^^^^ and F^^ 



(3d) 



Evaluating the singular pieces of -^"(3^,) ^^nd F^t^^-^ is simply performed by substituting in (6.4) and (6.5 1 and expanding 



to find 



+ 2(^A^(aM +pM-V,) j dr' Dr{z^^,z^'')Ir{z^') 

+ ^-a^lliz'^) + (a'' + P^'^V,) / dr' Dr{z^,z^^')Ir(z^^') 



(6.41) 



VII. RENORMALIZATION 



Combining our results for the first, second and third order self force corrections from (6.7l, (6.201, (6.26), (6.38) 



(6.40) and (6.41 1 gives the regularized self force through third order in e, 



A 



2m2, V47r 



2m4j \47rJ 



+ 



z\ 2m?c\c2 / A \ 3 



in J 



+ 



fca^AV 
47r I 



in 



m^c^c^ ( A 



m c\C2 3m C1C2 / A \ ra CiC^ ( A 



4 2 2 / /* /* 

+ ^^"^V.) + / c^r' ^?;^(^^ + 0{e^) . (7.1) 



The appearance of higher-order divergences (e.g., ^ A'^a^) and history-dependent divergences (e.g., ^ AP^'^W i,l'^{z^)) 
may seem unsettling. However, despite this seemingly pathological structure all the divergences through third order 
in e can be absorbed into purely local and time-independent counter terms in the effective action while the remaining 
terms independent of A cannot. 



More specifically, we add to the action in (5.14) the following counter terms 

S2 



dr 



mpi 



2mt 



dT^^{z^') 



(7.2) 



The Feynman rules for the counter terms above are shown in Figure [5] and their contributions to the effective action 
(and hence the self force) are given by the Feynman diagrams in Figure |6] Following similar steps we used in deriving 
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Cl (52 Ci (5^ C2 Ci 

FIG. 6: The counter term diagrams that contribute to renormalizing the self force through third order in e. 

and regularizing the self force through third order above one can show that the contributions of the counter terms to 
the self force are 



CT 



- S„ 



mciSi 



An J 



V.)/i?(z^)| 
1 



f 2mciSi 



An 



m^ciC2(5i 



-V 

An J 



A 

An 



Am?ciC25i 



'"-pi / "•'pi 

dr' Dn{z^',z^^')lR{z^')]\ - 



JL 
An J 



A 

An 

51 



2r7i^ciC2(5i 

^4 



(7.3) 



Notice that the divergent terms automatically organize themselves to have the same structure as the divergent terms 



appearing in (7.1 1. That is, the divergences appearing in both (7.1 1 and (7.3l are proportional to a'^, (a'^ +P^''V^)/ij 
and (a^ + P^'''V^)(/^/2 + •••)• This is no accident. In order for the divergences to be removed, the structure of the 
counter terms must exactly match the divergent parts of the self force. This can be viewed as a consistency check 
that our calculations to this point are correct. In addition, the terms in (7.1 ) not involving A cannot be absorbed by 
these counter terms. 

The divergences in _F^(r) are cancelled by requiring Sm, 6i and S2 to eliminate all terms proportional to A, which 
implies 



2^ 



A \ m'^clc2 
^ 2mt, 



A 

An 



Si 



S2 = 



nri, \An J 



-pi 

2^2 



771 c: 



A 

An 



"pi 



4- ? 9 
TO CJC2 

m CJC3 



-V 

An J 

A 
An 



A S 

m C1C3, 
6to2, 



-V 

An J 



m 



pi 



A 

An 



(7.4) 
(7.5) 
(7.6) 



Notice that there are three divergent terms appearing in the self force in (7.1 1 and there are three counter terms that 



absorb these divergences by renormalizing m, ci and C2. Had we computed the self force only up through second 



order in e then (7.1 1 would only have two divergent terms, proportional to and [a^ + P^'^V u)Ir, so that only two 
counter terms would need to be introduced {5m and 5i). The key point is that at ©(e") the self force will have n 
divergent terms that can be removed by introducing n counter terms into the action so that the theory can be made 
finite at any order. Therefore, the nonlinear scalar model is renormalizahle at any given order of e. This statement is 
independent of which field variables are used (0 or iji), since the actions are the same in either case. 

With all divergences cancelled by the appropriate counter terms, the renormalized self force through third order in 
e is given by 

F^(r) = (a^ + P^'^V.)(^/H(2^)- ^^^f;^4(z^)+ / dr' Dn{z^ , z^^' )In{z^^' ) 
i^pi ^pi V2 J 

^"^Pf/fl(^^) / dr' Dr{z'',z^'')Ir{z^')+ [ dr'dr" Dn{z'',z'^')DR{z''' ,z''")Ir{z^'") 



m 



m^CiCs ( 1 



/ dr'DR{z^,z^')ll{z^')^ +0(e')} , 



(7.7) 
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which is the main result of this paper. 

As stated before, in dimensional regularization all power-divergent integrals vanish so that A = 0. Therefore, 
whenever a regularized integral gives a term proportional to A we can simply drop that term from the computation. 
This is not the case with log-divergent integrals, which describe physically relevant screening effects. Instead, one must 
deal with these divergences as we have here by identifying the counter terms that cancel them. However, logarithmic 
divergences depend on the renormalization scale and induce a classical renormalization group (RG) running of the 
appropriate coupling constant. The possibility for classical RG running was first recognized in Ref. [30'. 

In the nonlinear scalar model presented here, all interactions occur between the field and the particle only (i.e., there 
are no field self-interactions) . Hence, the only kinds of divergences that can be generated are proportional to A^ for 
some integer p > 0, which vanishes in dimensional regularization. In other words, the theory described by the action 



in (5.141, or equivalently in (5.5), does not generate any logarithmically divergent integrals, only power-divergent ones. 



Therefore, in this model we can simply make the replacement 

D,,t{z^,7y)^ Dr{z'^,z'^') (7.8) 

wherever we see a retarded propagator in an expression connecting two points on the worldline ~ the resulting 
expression is regular. 

Dimensional regularization is certainly a convenient tool for regularizing the self force in the nonlinear scalar model. 
However, its use is necessary for gauge theories, e.g., for the gravitational self force. The reason for this is simply 
that the counter terms induced by the divergences are not necessarily restricted to be gauge invariant under small 
coordinate transformations on the background spacetime. For example, in calculating the finite-size corrections to the 
radiation-reaction force on an electrically charged extended body one must include counter terms such as / dr A^j^A^^ 
and J dra^A^, among others, which are obviously not gauge invariant [53J. Thus, using a cutoff regularization 
prescription, for example, would render the resulting theory finite but gauge violating. In dimensional regularization, 
the counter terms are zero for power-divergent integrals and the renormalized theory is finite and remains gauge 
invariant. The fact that dimensional regularization preserves the symmetries of the theory [70J is then quite crucial 
for computing the regular part of higher-order gravitational self force corrections. 



We could have obtained (7.7 1 with much less work had we simply used the replacement rule (7.8) in the formally 



divergent self force expression in (7.1). However, it is important to highlight the consistency of our renormalization 
program, especially since one may worry about the consistency of a point particle treatment for calculating higher- 
order self force corrections and/or about the subtleties associated with renormalizing gauge theories at higher orders, 
as discussed above in the previous paragraph. 

VIII. EQUATIONS OF MOTION THROUGH THIRD ORDER 

The regular self force equations of motion are given by ma'^ = where is given through third order in e by 



(7.7). All of the terms proportional to a'^ in this equation can be combined into one function rDw(^^) so that the 



effective mass of the compact object is mcff (r) = jtiFd-vvI-^^'') where 

„2 ^2„2 



izn-l~"4lRi^n + "^(llU^n+ f dr'Dniz^,z^')Iniz^^') 

/fl,(z^) J dT'Dniz^,z^')Iniz''' ) + J dr' dr" Dr{z^ ,z^^^ )Dr{z^\z^'" )Ir{z'^" ) ^ 

/ dr'DR{z^,z>^')ll{z^')^+0{e') . (8.1) 



m'^C{C3 ( 1 



2m6, 

In doing so, we find that the equations of motion can be written as 

mrDw(^^)a^ = -mP^''V^rDw(2^) (8.2) 
or, equivalently, dividing both sides by mFDw 

^ -PA'>'V^lnrDw(2'') • (8.3) 
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rupi 



3 1 suggests that FdwI-^^) = C'(V'-R/"^p^) when comparing with (5.201. If this is true, then the regular 



The form of 

part of ip in the Detweiler- Whiting scheme is inferred to be 

iPbXx) mci 



dTDR{x,Z^)lR{z^) 



m C1C2 



'i 9 



dTDR{x,z^)ll{z^) + Oie'') 



(8.4) 



In the next paper in this series we will show that (8.4 1 is indeed the correct expression for the i?-part of ip through 
0{e^) by explicitly calculating the radiative field emitted by the SCO through third orde r [29] . We will also compute 
the self force through third order from these emitted waves and find agreement with (8.3 1, thereby establishing a 
complete formalism that can self-consistently compute the SCO equations of motion and waveform. 

One can interpret mrDw(2^^) as an effective mass as in (8.2 1 or equivalently as an alteration to the self force 
experienced by the SCO as in ( |8.3| . If the effective mass of the SCO, mFDWi decreases under a constant self force, 
then the body is accelerated to higher velocities because there is effectively less inertia to resist the influence of the 
force. On the other hand, one could interpret the same dynamics as arising from a body with constant inertia but with 
a stronger self force acting on the body compared to the former case, which likewise makes the body move faster. The 
points of view (and hence where we place FdwI-s^^) in the equations of motion) give rise to the same dynamics despite 
the ambiguity in the interpretation. It is likely that such an ambiguity is also manifest in higher order gravitational 
self force corrections and probably originates from trying to localize the mass-energy of the SCO interacting with the 
background spacetime. 

Instead of splitting the retarded propagator into its regular and singular parts we could have equally chosen to 
decompose D,:ct{x,x') into its direct and tail parts as in Section IV We will not reproduce the steps for deriving and 
renormalizing the self force through third order in the latter decomposition since they are similar to the manipulations 
presented in Sections |IV| and |VI[ After a straightforward but tedious calculation we find 



meff(T)a^ = P^'' 



mz 



9 



4 2 2 

m CJC2 



Uiz^^) + Il-'\z^) 

2Mznit.-A{zn + ii^'\znit.ii{zn + nm 



0+ 



dT'V,i5rct(;2^,^^')/tail(2^ 



"pi 



J/tail(z^) 1 tail,,M^ -P^tail(^ 

dT An " ^ ' dT 



+ /*-'(z^) lim / dr' D,,,{z^,zy)I,,n{z^') + I,,iiizn lii-nl dr' V,7^,et(^^ ^^')/taii(2''') 



0+ 



2/,(z'^) lim 



dT'D,,^{z^^,z^'')I^,■a{z^') 



— lim 
47r e^0+ J^r 



dT'^,D,,,{z>',zy) 



dr' 



2ml, 



lim / dr"^ ^.D.Az^" , z^' ) lim 



^f.lLA^n + C^^iznilsii^n + lim 

e-i.O+ 



dr" D,,t{z^ ,z^ )/taii(^^ ) 
' ' V.A■et(^^;^''')/tlil(^''') 



+ 0(e4) 



where Itai\{z'^) and /*^''(z^) are the tail integrals 

/tail(z'^)= lim / dT'Z?rct(2^2^') 

/^^"(z'^)= lim / dT'W,D,,t{z^,z^'), 
which are defined only on the worldline, and 



(8.5) 

(8.6) 
(8.7) 
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is the local radiation reaction force in a non- vacuum background spacetime. The effective mass of the compact object 
is meff(T) = mrHad(^'^) where rHad(2;'') is given by 



rHad(z^ 



4 9 

m c^C2 



S 9 



Il^izn + hm 

£-i.0 + 



/taii(z'') lim 



dr' D,,t{z^,z^ )/taii(^M 
— lim / 

47r e^0+ 



lim 



47r e^0+ J_oo 

dT'Act(^^z^') lim / dr"Act(^^',^^")/taii(^^") 



'.^/tail(^^ ) 



dr' 



m^cf C3 
2mL 



tail 



(z^) 



lim 



(8.9) 



Writing the self force in terms of the tail part of the propagator does not provide a condensed form as in (8.3 1 since 



the right side of (8.5) is not proportional to Vj^rHadl-z^)- This is because of the appearance of /^(z'^), which is not 
proportional to Vjy, and because the covariant derivative can not be pulled through the tail integrals in (8.51 without 



generating other terms that are not proportional to V^,. However, the form in (8.5 1 and (8.9 1 clearly shows the role 



of tails-of-tails and tails-of-tails-of-tails in the self force and effective mass, respectively. 



The self force in (8.5 1 also depends explicitly on Da^/dr so that an order reduction procedure should be implemented 



to exclude unphysical runaway and acausal solutions from appearing in the theory - order reduction is a standard 
treatment of this issue in theories describing radiation reaction particle dynamics (see e.g., |71|). 



IX. CONCLUSION 



In this paper we introduced a nonlinear scalar model for extreme mass ratio inspirals that is a natural analogue 
of the corresponding perturbative General Relativistic description. This model should be useful for studying the 
role of higher-order self force corrections for building high-accuracy waveforms for precise parameter estimation, for 
quantifying the effect of transient resonances on the phase evolution of waveforms, and for providing a sufficiently 
simple context to develop numerical methods for computations of higher-order self force corrections that may then 
be applied to the gravitational problem. This last program should be particularly useful for calibrating semi-analytic 
models [17J, for building phenomenological-based hybrid waveforms [18J, and for constructing gauge-invariant EMRI 
observables at higher orders in e. 

This model has the interesting feature that, despite being constructed to be as similar as possible to perturbed 



General Relativity, one can always perform a field redefinition using (5.10) so that the initially nonlinear scalar field 
theory is transformed to a linear one (with suitable changes in the field-particle interactions) . Such a transformation 
provides a very clean and simple derivation of the self-force corrections at higher orders since all contributions coming 
from the self-interaction of the field are subsequently removed. It is natural then to wonder if such a transformation 
can be made when describing gravitational EMRIs to remove some, if not all, of the nonlinear self-interactions of 
the metric perturbations. Indeed, it may be that a combination of a field redefinition and a gauge transformation 
(one that is O(e^) so as to keep the linearized metric perturbations in the Lorenz gauge at leading order) successfully 
removes all nonlinear self-interaction terms at a given order. If this is the case then calculating higher order self force 
corrections will be much easier to derive analytically and solve numerically. We will discuss these issues further in a 
future paper. 

In order to calculate the self force on the SCO we demonstrated how the usual action principle fails to retain 
the outgoing boundary conditions when integrating out the field - a conservative dynamics for the worldline ensues 
because the full Lagrangian is time-reversal invariant. However, using the difference in the actions for two worldline 
and field histories (one evolving forward in time and the other backward) allows one to "break" the time-reversal 
symmetry naturally contained in the usual action principle and incorporate the time-asymmetric outgoing boundary 
conditions into the effective action. In this way, one can describe the open classical system dynamics of the SCO from 
an action principle and derive the self force (and waveforms - see the next paper in this series) in a self-consistent 
manner. 

Using this new action principle for open classical systems together with the effective field theory formalism we 
derived the self force on the SCO through third order in e. Our main results are given in (7.7) and (8.3). We 



separated the singular parts of the self force from those that are regular on the worldline using the Detweiler- Whiting 
decomposition for the retarded Green's function. For completeness and comparison, we also presented the more 
complicated third-order self force expressions in terms of Hadamard's tail propagator in (8.5) and (8.9). Despite the 
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FIG. 7: The relation between the worldline (solid line) and a geodesic (dashed line) intersecting the worldline at the coordinates 
2''(r) and z'^ir') for constructing the quasilocal expansions of various quantities. 



appearance of higher order poles and history-dependent divergent forces - see (7.11 - we found that only three local and 
time-independent counter terms are needed to absorb all of these divergences, which amounts to a renormalization 
of the SCO mass m and the coupling constants ci and C2 appearing in (5.14). In addition, we showed that this 
nonlinear scalar model only generates power-divergent integrals that automatically vanish in a regulariziation scheme 
that manifestly preserves the symmetries of the theory (e.g., dimensional regularization) . No logarithmic divergences, 
which represent physical effects of screening by the scalar perturbations, can be generated in this model and thus 
none of the coupling constants appearing in (5.141 or, equivalently, (5.51 undergo a classical renormalization group 
running. 

Our results in (8.3) suggest an ambiguity in defining the effective mass of the SCO since, in this model, the 
effective mass is proportional to the SCO's inertial mass (moff = Fdw™) ^^nd one can associate the proportionality 
factor FdwI-z'') with a contribution to the effective mass or with an additional contribution to the self force. Both 
interpretations yield the same acceleration on the SCO and the same waveform and physical observables. 

In the next paper in this series we shall derive the expressions for the scalar perturbations emitted by the SCO 
through third order in e. We shall also show that one can more easily derive the self force by calculating the regular 
part of the radiative field and then evaluating it in the worldline equations of motion (see the discussion below (8.3)). 
This suggests that current 3+1 self force codes [22 may need only minimal modifications to compute higher-order 
self force corrections. In regularizing the third-order waveforms (which also contain divergences since the field is 
sourced by an effectively point-like object) we will show that exactly the same counter terms used in regularizing 
the third-order self force expressions - see (7.4l-(7.6l - also regularize the scalar perturbations and the self force 
expressions derived from them. Therefore, the use of point particle treatments for the SCO in conjunction with our 
renormalization program will be explicitly self-consistent. 
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Appendix A: Quasilocal expansions 

In this Appendix we provide the quasilocal expansions of various quantities used to regularize the expressions for 
the self force on the particle. 

For two points z'^ and z'^ on the worldline the quasilocal expansions of the direct /singular part of the retarded 
propagator ~ A^/^S{cr) follow from the expansions near coincidence (when s = r' — r is zero) of the van Vleck 
determinant A{x,x') and Synge's world function a{x,x'). The quasilocal expansion of the covariant derivative of 
Synge's world function, cr,y = cr-i, can be derived by noting that a unique geodesic (with coordinates a;''(A) and affine 
parameter A) intersects the accelerated worldline at the points z^(t) = x'^{0) and z''(t') = x^{X); see Figure [7| For 
r' near r and A near it follows that a;^(0) — a;^(A) = z^(t) — z^{t') implies 

-Aa;'''(0) - ya;"^(0) + • • • = -si^(r) - y^^M + • • • (Al) 
where a ' or ' denotes differentiation with respect to A or r, respectively. The geodesic equation can be used to 
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eliminate terms with two or more parameter derivatives since 

x"^{X) = -K,{x''{\))x"^{X)x'^{X) 



(A2) 



Finally, using the fact that is proportional to the tangent on the geodesic at A = so that au{x^,x^ ) = — Aa;J^(0) 
[57| gives 



Da 



2! " 3! dT 

Sygne's world function is constructed from the identity 2(t = a^a^ , 



(A3) 



and the quasilocal expansion of (5(ct) then immediately follows, 

8{s) 



\da jds 
S{s) 
\s\ 



(A4) 

(A5) 
(A6) 



The covariant derivative of 6{(j) in the quasilocal expans ion follow s by noting that ~ (J^{ds/da)d/ds. One can 
show that ds/da ~ (da/ds)^^ and use the expressions in (A3 1 and (A4| to find that 



V.5(a(z^z''')) = a,^-^<5(a(z^z^')) 
da ds 

s s2 23a" 
2 dr 



+ 0(.^))^<5(a(z^z^')) 



(A7) 
(A8) 



The quasilocal expansions of the van Vleck determinant and its covariant derivative follow from its expansion near 
coincidence [72J, 



12 
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Inserting (A3 1 into these expressions gives their quasilocal expansions. 



6 



(A9) 
(AlO) 

(All) 
(A12) 



These are all the expansions used in regularizing the self force corrections in this paper. 



Appendix B: Dimensional regularization 

We have mentioned several times in the main body that the divergent integral A vanishes in dimensional regular- 
ization. We will prove this here. 
First, note that 



and from (A6l it follows that 



A=i/ dsS{aiz'^,z^-)) . (B2) 
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FIG. 8: The Feynman rule for the interaction of n scalar fields on the worldline. 



[In fact, it should be clear from (6.4 1 that A = 47r J 
-s^/2 where s 



oo 

oo 



ds Ds{z^ , ).] In Fermi normal coordinates, a{ 



r' — T. The integral only has support when s ^ so that we can introduce a tangent space at the 
point z^{t) to define a Fourier representation for the integral in (B2|. In particular, if k'^ is a "momentum" vector in 
this tangent space then 



A = i hm 

2 nij,— i-O 



ds 



d'^k 



2 ™™oy_^ (27r)3k2 



^ (27r)4 k^ + ml 
d^k 1 



ra. 



(B3) 
(B4) 



where we have included a small mass in the denominator to make the integral well-defined when k^ — ^ 0. 

The divergence in evaluating the integral (B4 1 comes from those momenta with arbitrarily large magnitudes k — ^ oo. 



For a small enough spacetime dimension the ultraviolet behavior of A can be improved so that the integral is finite. 
This is the motivation for dimensional regularization wherein the spacetime dimension is parameterized by a complex 
number d so as to find a regime in the complex d-plane where the integral is finite. After evaluating the integral there, 
the result is analytically continued back to the physical four spacetime dimensions. 
Here, the divergent integral (B4) in arbitrary dimensions is 



A = 



1 



4-d 



2 {2t:Y' 



dfl lim 

d-2 mj,->0 



dk 



k 



d-2 



rrij 



where /i is an arbitrary parameter that preserves the units of A. Focusing on the last integral we see that 



lim 



dk 



fc2 



lim 



'4,-^0 2 ' 



(B5) 



(B6) 



which obviously vanishes in the physical four dimensional spacetime. Therefore, A = in dimensional regularization 
implying that one may replace D^ct{z^ , z^^ ) by Dfj{z^, z^ ) to obtain the correct regular self force, waveform, etc. as 
discussed in Section VII near (7.8 1. 



Appendix C: Feynman rules 

The effective field theory approach [30_ utilizes Feynman diagram techniques that are useful for efficiently obtaining 
the effective action at a given order in perturbation theory. Specifically, writing down the connected and tree-level 



diagrams at a given order (see Section V B I is equivalent to solving the wave equation through that order and 
substituting that solution back into the action to get the contribution to the effective action at that order. 

There are rules, called Feynman rules, for translating Feynman diagrams into mathematical expressions. The 
Feynman rules needed for computing the self force, or any other quantity in the theory, are given by: 

• For each worldline vertex represented by a dark circle and labeled by ai , • • • , a„ and c„ write down a factor of 
Tai---a„{x] z^]. These are the worldline vertices. 

• For each curly line write down a factor of D°'''{x,x') connecting worldline vertices labeled by the indices a and 
b at spacetime points x and x' on the worldline. 

• Sum over all labels oi • • • a„ in the chosen representation (e.g., the ± basis) and integrate over all spacetime 
points. 



• Divide by the appropriate symmetry factor of the diagram. 
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The symmetry factor for a given diagram is found by counting the number of ways one can permute the histories' 
indices (i.e., ai, 02, . . .) while retaining the same set of propagators appearing in the diagram. The vertex Taj^...a„{x; z] 
comes from the interaction terms in the full action (5.141 involving a worldline integral over n scalar fields, 



1 



where 



T{x;z'' 



mc„ 



m 



dTTp^'izi') = — / ^^"(a;)r(a;;z^ 



pi 



3I/2 



(CI) 



(C2) 



-pi 



which also defines V{x; z'^]. The indices ai, • • • , a„ on Ta-^...a,X^) A come from doubling the dynamical variables to 
ensure that the variational principle and the effective action are compatible with the outgoing boundary conditions 
on the scalar field as discussed in Section [III B[ Specifically, the contribution of this interaction term to that action is 



1 



i^r(a;)T(a;;zi] 



i^'i{x)T{x-Z2] 



,p^^{x) ■ ■ ■^^"{x)ca,...a„''Tb{x;z] 



where Tb{x; z] = T{x; zb\, capital Roman letters take values in {1, 2} and 



1 Ai = ...= 

otherwise 



Ap = 1 
Ap = 2 



(C3) 



(C4) 



In the ± basis, Ca^-.-aJ' = ^ai^^ ■ ■ ■ -^a^^" bcAi---a,^^ where Aa^ is given in Section IIIC 

Ta^...a„ix;z] = -''^^^Ca^...aJ'Vb{x;z\ 



Therefore, 



(C5) 



where Vf,(a;;z] = Ki,^V{x;zb\- This worldline vertex is represented by the diagram in Figure [s] and is all that is 
needed to construct the self force through third order in e in this nonlinear scalar model. 
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